A model of deterministic detector with dynamical decoherence 
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We discuss a deterministic model of detector coupled to a two- level system (a qubit). The detector 
is a quasi-classical object whose dynamics is described by the kicked rotator Hamiltonian. We show 
that in the regime of quantum chaos the detector acts as a chaotic bath and induces decoherence of 
the qubit. We discuss the dephasing and relaxation rates and demonstrate that several features of 
Ohmic baths can be reproduced by our fully deterministic model. Moreover, we show that, for strong 
enough qubit-detector coupling, the dephasing rate is given by the rate of exponential instability of 
the detector's dynamics, that is, by the Lyapunov exponent of classical motion. Finally, we discuss 
the measurement in the regimes of strong and weak qubit-detector coupling. For the case of strong 
coupling the detector performs a measurement of the up/down state of the qubit. In the case of 
weak coupling, due to chaos, the dynamical evolution of the detector is strongly sensitive to the 
state of the qubit. However, in this case it is unclear how to extract a signal from any measurement 
with a coarse-graining in the phase space on a size much larger than the Planck cell. 

PACS numbers: 05.45.Mt,03.65.Ta,03.65.Yz 



I. INTRODUCTION 

The long-standing problem of quantum measurement 
Qi H> Hi has recently gained a renewed interest due to its 
relevance for quantum information processing. Indeed, 
one of the requirements for the physical implementation 
of quantum computation is the ability to readout a sin- 
gle two- level quantum system (qubit). This problem has 
been solved in ion-trap quantum computation '3,'5j using 
quantum jump detection. In solid-state implementations, 
the single-qubit measurement is very chal leng in g and has 
been widely discussed H S H H Efl [H El [13. More- 
over, various readout schemes have been experimentally 
realized [11 El E [13 . 

A detector can be seen as a complex quasi-classical 
object coupled to a quantum system. It is therefore in- 
teresting to investigate the dynamical evolution of con- 
crete system-detector models. In this paper, we intro- 
duce and study a model in which the detector exhibits 
a complex quasi-classical dynamics and is coupled to a 
two-level quantum system. The detector is described by 
the kicked rotator Hamiltonian 0|, so that its dynam- 
ics can be, in the classical limit, chaotic, integrable or 
with mixed phase space. An interesting feature of our 
model is that it can be, in principle, realized by means 
of cold atoms in optical lattices. Cold atoms exposed to 
periodic standing waves of light proved to be an ideal 
testing ground to explore the quantum dynamics of non- 
linear systems ^19^, i2Qj and, in particular, made possi- 
ble to realize experimentally the quantum kicked rotator 
[2ll[2l[2l[2l|. Moreover, the sensitivity of the quantum 
dynamics to a control parameter has been recently mea- 
sured for a classically chaotic system and suggested as a 
promising technique for precision measurements [25| . 



It is clear that the process of quantum measurement 
involves a tradeoff between information gain and distur- 
bance ^2^. This means that the price to pay for the 
information to flow from the system to the detector is 
the backaction by the detector. The backaction disrupts 
the coherent evolution of the quantum system. This en- 
tails a fundamental connection between decoherence and 
quantum measurements. In this context, the detector 
model introduced in this paper can be seen as a model 
of chaotic environment |23|. We will show that indeed 
the coupling to the detector induces decoherence in the 
quantum system. We will also show that the dynamical 
decoherence resulting from our fully deterministic model 
has properties very similar to the decoherence induced by 
a suitable, dissipative coupling to a non-deterministic en- 
vironment. An interesting trait of our model is that the 
decoherence rate is given, for strong enough coupling, by 
the Lyapunov exponent, namely by the rate of exponen- 
tial instability of the detector's dynamics in the classical 
limit. 

An interesting problem is what is the system-detector 
coupling strength required to measure a single qubit sys- 
tem. We will show that, in our model, a strong coupling 
allows us to measure the up/down state of the qubit. On 
the other hand, due to chaotic dynamics, also a weak 
coupling induces significantly different dynamical evolu- 
tions of the detector in the case of qubit up or down state. 
However, it is not clear how to extract a detectable signal 
from this difference. 

The paper is organized as follows. In Sec. II, we in- 
troduce our deterministic detector model and describe 
its possible implementation by means of cold atoms in 
laser fields. In Sec. Ill, we discuss the decoherence of 
the dynamical system, induced by the the coupling to 
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the detector. The results are compared with those of 
a phase damping map for the system's density matrix, 
describing the decoherence process in the quantum op- 
erations formaUsm. In Sec. IV we analyze the depen- 
dence of the detector's response on the system-detector 
coupling strength. Our conclusions are given in Sec. V. 
We also present an alternative derivation of the phase 
damping map, based on the master equation approach 
(Appendix A) and discuss the continuous limit of this 
map (Appendix B). 



II. DETERMINISTIC DETECTOR MODEL 

We consider the interaction of a single qubit with a 
quantum kicked rotator. The overall Hamiltonian H 
reads as follows: 

H = Hs + Hd + Hint , 
Hs = Sa-j;, 

Hd^\ + Kcos{e)Y,5{t~m), (1) 

m 

Hint = ecCTz cos(^) (5(t - m). 



Here Hs denotes the single-qubit Hamiltonian, Hd the de- 
tector's Hamiltonian (kicked rotator model 1^) and Hint 
the qubit-detector coupling. The Hamiltonian Hs in- 
duces Rabi oscillations with frequency uji = 25 between 
the qubit levels |0) and |1) (|0), |1) denote the eigenstates 
oicFz corresponding to the eigenvalues -f 1 and —1, respec- 
tively). The detector is a particle moving in a periodic 
potential switched on/off instantaneously (kicks) at time 
intervals At = 1. The time t in Eq. is measured in 
number of kicks. We have [p, 6] = —ih, where h is the 
effective dimensionless Planck constant. The prop erties 
of the quantum kicked rotator are described in . The 
interaction Hamiltonian is also kicked with the same time 
period At — 1. Wc note that the interaction Hamilto- 
nian (cx ffz) does not commute with the system's Hamil- 
tonian (cx (Ta;), namely we do not discuss nondemolition 
measurements 26]. The unitary operator describing the 
evolution of the overall system (qubit plus detector) in 
one kick is given by 



U = exp 



K 



■ COS( 



exp 



^2h 



exp(-i(5(72;). 



(2) 

Therefore, the effective kicking strength Kcs = K + ecCz 
depends on the up or down state of the qubit. For this 
reason, the kicked rotator can be, in principle, used as a 
detector. We define e = ec/h as the coupling strength in 
units of h. 

The classical limit for the detector corresponds to 
h 0, while keeping K constant. The classical dynamics 
is integrable at K = and exhibits a transition to chaos 
of the KAM (Kolmogorov-Arnold-Moser) type when in- 
creasing K. The last invariant KAM torus is broken for 



K \ and the motion is full y ch aotic (with no visible 
stability islands) when K ^ 1 |23|- In this paper, we will 
consider the detector as a quasi-classical object, that is, 
we will take ft ^ 1. 

In the numerical simulations discussed in this paper, 
we consider initial separable states l^*) = IV's) ® 
where \tps) = +/3|1) is a generic single-qubit state 
and \4)d) is a Gaussian wave packet describing the initial 
state of the detector (centered at p = 0, = vr, with 
the uncertainties along 9 and p proportional to Vft) • We 
consider the kicked rotator dynamics on the torus < 
9 < 2tt, —it < p < tt. The number Nd of quantum levels 
describing the detector's dynamics is Nd = 27r/ft. We 
consider 2"^ < Nd < 2^^, corresponding to 4.91 x 10^^ > 
ft > 1.92 X 10"^. 

To close this section, we would like to comment on 
the possibility to implement, at least in principle, our 
model using cold atoms in a pulsed optical lattice cre- 
ated by laser fields. For this purpose, the atoms should 
be prepared in a given hyperfine level. A microwave ra- 
diation then creates Rabi oscillations between this level 
and another hyperfine level 24]. These oscillations are 
described by the single-qubit Hamiltonian Hs- Hence, 
the two hyperfine levels correspond to spin up and spin 
down. The kicking strength K depends on the hyper- 
fine level, as required in our model. After a sequence of 
kicks, the lattice is switched off and the momentum dis- 
tribution of the cold particles can be measured. As we 
will show in this paper, this can give information on the 
state (hyperfine level) of the qubit. We point out that 
a similar experiment has been proposed |29j and experi- 
mentally realized 124"] with a different purpose, namely, to 
measure the fidelity of quantum motion under a Hamilto- 
nian's perturbation (see, e.g., [s^, |^ IH, 113, 0, 113 and 
references therein). Therefore, the experimental imple- 
mentation of our detector model is in principle possible, 
even though the realization of the kicked rotator dynam- 
ics has been so far possible only for values ft 1 or larger 



III. DYNAMICAL DECOHERENCE 

In this section, we study the dynamical decoherence of 
the system, induced by the couphng to the detector. In 
this context, the detector can be seen as a microscopic 
model of a chaotic bath. We consider two time scales Ti 
and T2 that characterize the relaxation of the diagonal 
elements of the system density matrix and the decay of 
the off-diagonal elements, respectively. As we are inter- 
ested in the measurement of the spin polarization along 
the z axis, the time scales Ti and T2 will be derived in 
the basis of the eigenvectors of CTz. A short discussion of 
the time scales that characterize the evolution of diagonal 
and off-diagonal elements of the qubit 's density matrix in 
the basis of the energy eigenstates of the qubit (the ax 
basis) will be given in Appendix B. We consider the case 
(5 <C 1, corresponding to a time between kicks (qubit- 
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FIG. 1: Dephasing of the spin system for K = 8, e = e^/h — 
0.3, 5 = 0.2, h = 7.67 x 10"", and initial spin state |V's) = 
(10)+2|l))/v^. The straight line fit gives I Poll =aexp(-r2t), 
with a = 0.40 and r2 = 2. 17 x 10"^. Here and in the following 
figures the logarithms are decimal. 

detector interactions) much smaller than the period of 
the free-qubit Rabi oscillations '36]. We also derive, in 
the framework of the quantum operations formalism, a 
phase damping map for the system and compare the re- 
sults of this approach with those obtained from the exact 
numerical solution of the overall dynamics (system plus 
detector). 

The dephasing induced by the detector is illustrated 
in Fig. n We can clearly see the exponential decay of 
the non diagonal term poi of the reduced density ma- 
trix p{t) = Tr<j[|*(t))(^'(t)|] of the two-level system (Tr^ 
denotes the trace over the detector). Since the overall 
Hilbert space is finite, the exponential decay is possible 
only up to a finite time, after which quantum fluctua- 
tions determine the residual value of poi- Nevertheless, 
the decay rate r2 can be clearly extracted from a fit of 
the short-time decay 

IpoiWI «exp(-r20. (3) 

This gives us the dephasing time scale T2 — 1/T2- We 
also note that the exponential decay is superimposed to 
Rabi oscillations with period given by tn = 2tt/ll}ji — 
tt/6. This is due to the fact that the Hamiltonian Hg 
induces free rotations of the spin around the a;-axis with 
the rotation period being tn. 

The dependence of the decay rate r2 on the coupling 
strength e ~ ec/fi is given in Fig.|3 for different values of 
the effective Planck constant h. For weak coupling e ^ 1, 
we have in average 

T2 ~ eV2, (4) 

in agreement with the expectations of the Fermi golden 
rule. For e > 1, the dephasing rate saturates to an e- 
independent value. The discussion of the strong coupling 
regime is postponed to the end of this section. 

As we have seen in Fig. ^ at long times poi{t) oscil- 
lates around a residual value pros- In Fig- 3, we show 
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FIG. 2: Dependence of dephasing rate r2 on for S = 0.1, 
h = 4.91 X 10"^ (plus), ;i = 1.23 X 10"^ (triangles), h = 
3.07 X 10"^ (circles), and h = 7.67 x 10"* (stars). The initial 
state and K are the same as in Fig.Q The straight line shows 
the average fit r2 — Ae^, with A = 0.57. 

the dependence of prcs on the overall size A'^ = NdNg of 
the Hilbert space, where Nd and Ng = 2 are the dimen- 
sions of the Hilbert spaces for the detector and the spin, 
respectively. We can see that 

Pres oc 1/Vn, (5) 

as expected from the following simple statistical estima- 
tion. We assume that at long times the overall wave 
function is ergodic, that is, 

1 Na 

i*w) = EE^»wi*)i")' (6) 

i=0 n=l 

where the coefficients Cm(t) have amplitudes ^ 1/\/A 
(to assure the wave function normalization) and random 
phases. Therefore, poi{t) = X^n conC^^ ~ 1/Vn (sum of 
N/2 terms of amplitude ~ 1/N with random phases). 

The relaxation of the diagonal componen t pi 1 (t) to the 
asymptotic value pn = ^ is shown in Fig. ^[33. It can be 
seen that the population relaxation is exponential, with 
rate Fi. Indeed, the evolution in time of pii{t) — ^ is 
well fitted by the curve asm(bt + <j)) exp(— Fit), with h = 
0.404 m 25 = 0.4 frequency of the free-qubit oscillations 
and Fi ~ 4.36 x 10"'^. This fit allows us to extract the 
relaxation time scale Ti = 1/Fi — 22.9. 

The dependence of the relaxation rate Fi on the cou- 
pling strength e is shown in Fig. |S1 At small e, the be- 
havior of the relaxation rate is similar to that of the de- 
phasing rate: we have 

Fi « F2 « eV2. (7) 

Moreover, the relaxation rate exhibits an interesting non 
monotonous behavior: it has a maximum at Fi ~ (5 and 
then decays when increasing e. This phenomenon is a 
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FIG. 3: Dependence of pros (obtained after averaging |poi(t)| 
over the time interval 500 < t < 2000) on the size A'' of the 
overall Hilbert space. Parameter values are the same as in 
Fig. The straight line gives pres = 1/VN. 




FIG. 4: Population relaxation, that is, evolution in time of 
pii (solid curve). Parameter values are as in Fig. The 
dashed curve shows the fit pn = | + asin(6t + (f)) exp(— Fit), 
with a = 0.5, b = 0.404, 4> = 0.405 and Fi = 4.36 x IQ-^. 



manifestation of the quantum Zeno effect pM l39| : re- 
peated measurements performed by the environment (the 
detector) prevent the system from relaxing. In the Zeno 
regime, the stronger the dephasing, the slower the relax- 
ation. In our model, this is true in a rather broad range 
of e, that is VS ^ e < 1 (see Appendix B). In the case 
of Ohmic dissipation one expects [Till I4(TI | 



(8) 



This is in a good agreement with our numerical data, 
shown in Fig. In the same figure, we compare these 
data with the theoretical curve Fi = B/e^, with the fit- 
ting constant B = 2.7x 10~^ (note that here S'^ = 10"^). 
For strong interactions e > 1 the relaxation rate exhibits 
an oscillating behavior. 

It is interesting to compare the results obtained from 
our deterministic detector model with those of a quantum 




FIG. 5: Dependence of the relaxation rate Fi on the coupling 
strength e, at 5 = 0.1, iC = 8, ft = 4.91 x 10"^ (plus), h = 

-3 



(triangles), H = 3.07 x 10 ^ (circles), and H 
"* (stars). The initial state of the qubit is \ips) 



1.23 X 10" 
7.67 X 10" 

(10) + 2\l))/^/5. The straight lines show Fi = Ae^, with 
A = 0.56 (solid line) and Fi = B/e'^ with = 2.7 (dashed 

line) . 



map derived for the system density matrix in the frame 
of the quantum operations formalism. 

The evolution in one period of time of the system den- 
sity matrix can most conveniently be written in the Bloch 
sphere representation, in which the coordinates x, y, and 
z are related to the matrix elements of p as follows: 
pii = ^(1 — z) and poi — -^{x — iy). The free evolu- 
tion between two consecutive kicks is ruled by the free 
Hamiltonian ffg, that is, 



From this equation we obtain 



X ~ X, 

y ~ cos{2S)y — sin(2(5)z, 
z = sin(25)y + cos(2^)z. 



(9) 



(10) 



the coordinates i, y, z corresponding to the density ma- 
trix p. 

We model the effect of the interaction with the detector 
as a phase kick. That is to say, we assume that in the 
interaction Hamiltonian of Eq. iQJ the angle Q is drawn 
from a random uniform distribution in [0, 27r]. Therefore, 
the density matrix p, obtained after averaging over 0, is 
given by 



P = 



1 

2^ 



deR{9)pR\e), 



where 



R{B) 



ie cos 9 

e 





le cos 6 



(11) 



(12) 
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For e <C 1, we obtain 

R{e) « / - iecos{6)a, - — 0082(6*)/. (13) 

Using this approximation, we end up with the map 
p = p + (^azpaz < cos^ 9 > -f? < cos^ 9 > p 

where in the random phase approximation we take < 

008^6' >— J- d9 cos^ 9 — i. This map is known as 

I — 1 1 — ^ 

the phase damping noise channel |41I l42l|. 

Map (dH) can be written in the Bloch sphere coordi- 
nates as fohows: 

X = (1 — e2).T = (1 — e^)x, 

y = (1 - = (1 - e')[cos{25)y - sin(2^)z], (15) 
z = z = sin(2(5)j/ + cos(2(5)z, 

where the coordinates x, y, z correspond to p. An alterna- 
tive derivation of map H15|l . based on the master equation 
approach, is provided in Appendix A. 

This map gives the evolution of the Bloch sphere coor- 
dinates in one kick and can be iterated. From the values 
of the coordinates x, y, and z after t map steps we can 
obtain poi s-iid pn at time t. An example of numerical 
solution of map p5|l is shown in Fig. IHl The exponential 
decay of poi and the exponential relaxation of pn to the 
steady state value ^ can be clearly seen. Similarly to 
Fig. Gland Fig. 01 we also have oscillations with the Rabi 
frequency. We also checked that the relaxation rates Fi 
and F2 are both proportional to e^. More precisely, the 
numerical data for the model (|15|l give Fi w r2 ~ O.SSe^. 
These rates are in good agreement with those from the 
numerical data obtained for the whole system (qubit plus 
detector) and shown in Fig. El and Fig.O Therefore, it is 
noteworthy that our fully deterministic dynamical model 
can reproduce the main features of the single-qubit de- 
coherence due to a heat bath. Our results are also in 
agreement with the analytical solution obtained in the 
continuous time limit and described in Appendix B. 

To close this section, we discuss the case in which the 
coupling strength is classically strong, namely ec ~ 1, 
corresponding to e » 1 for a quasi-classical detector 
{H 1). As shown in Fig. [7| the dephasing rate is given 
by the Lyapunov exponent A of the underlying classical 
chaotic dynamics of the kicked rotator. We point out 
that, in this regime, the dephasing rate is independent of 
the interaction strength. This result can be understood 
as follows. Assuming that the period of the free oscil- 
lations of the system is much larger than the dephasing 
time, namely 1/S ^ X, we have 

PoiW ~ poi{o){Miui+Y{Ud^-Ym, (16) 

where Ud,+ and Ud.- are the one-kick evolution opera- 
tors for the detector when the effective kicking strengths 
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FIG. 6: Dephasing and relaxation in the phase damping map 
HSU . Main figure: time dependence of |poi| for e = 0.225 and 
other parameter values as in Fig. The straight line fit gives 
the exponential decay |poi| = aexp{—T2t), with a = 0.31 
and r2 = 2.77 x 10~^ = 0.55e^. Inset: same as in the main 
figure but for pn. The behavior of pn (dots) is very well 
reproduced by the fit pn = ^ -I- asm{bt + 4>) exp(— Fit), with 
a = 0.3, 6 = 0.4, (j) = 1.52 and Fi = 2.30 x 10"^ = 0.45e^ 
(solid curve). 



are iiToff = K+ — K + ec and K^s — K- = K — ec, re- 
spectively ^2 . For = K_ , the initial Gaussian wave 
packet is centered at a stable fixed point and therefore, 
for short times, {Ud,-Y\4'd) ~ I'^rf)- On the other hand, 
the same fixed point is unstable for K^ff = K+. There- 
fore, in the quasi-classical regime and for strong enough 
perturbations the wave packet spreads along the direction 
of instability for the classical motion, with rate given by 
the Lyapunov exponent. We note that this phenomenon 
is closely related to the Lyapunov decay of the fidelity of 
quantum motion in chaotic systems ^ll [S^, |33|, (sJi |35| . 
Indeed, Eq. (|16|l shows that poi{t) is proportional to the 
fidelity amplitude f{t) = {^dml+YiUd-Ym. This 
quantity measures the stability of quantum motion un- 
der perturbations. More precisely, f{t) is the overlap of 
two states which, starting from the same initial condi- 
tions, evolve under two slightly different Hamiltonians. 
The above discussion shows that the fidelity decay has 
an interesting interpretation in terms of dephasing of an 
appropriate two-level systems. 



IV. DETECTOR MODEL 

In this section, we study the efficiency of the detector in 
the regimes of weak and strong system-detector coupling. 

Let us first discuss the strong measurement case with ec 
large ec ~ 1. In this case, the effective coupling strength 
KeS = K + ecCTz significantly depends on the up or down 
state of the spin system. Therefore, it is easy to find 
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FIG. 7: Time dependence of |poi| for K = 4.5, Ec = 
e27r/2""^ = 0.8, and 5 = 0.01. The initial state of the qubit is 
\tps) = (|0) + 2|l))/\/5, the initial detector state is a Gaussian 
wave packet with area size h centered at p = 0, ^ = tt. Data 
are shown for h = 4.91 x 10"^ (plus), h = 1.23 x 10"^ (trian- 
gles), h = 3.07 X 10"^ (circles), h = 7.67 x 10"" (stars), and 
h = 1.92 X lO"** (diamonds). The straight line represents the 
exponential decay with rate given by the Lyapunov exponent 
A«ln(jr/2) =0.81 m. 



regimes in which the response of the detector clearly de- 
pends on the state of the system. An example is shown 
in Figs. I8I9I In these figures, we have K = 4.5, while 
Kcs = K — Ec = 3.8 when the rotator is coupled to a 
down spin (cr^ = —1) and K^s = K + £c = 5.2 when 
the rotator is coupled to an up spin (cr^ = +1). The ini- 
tial Gaussian packet is centered at the fixed point p — 0, 
9 ~ TT. A linear stability analysis of the classical dy- 
namics of the kicked rotator shows that this fixed point 
is stable for < K < 4 and unstable for K > 4 and 
K < 0. Therefore, in the case of coupling to a down 
spin, the fixed point is stable {Kcs < 4), while the same 
point is unstable when the detector is coupled to an up 
spin {KcB > 4). In Fig.|Sl we see that, if the qubit is in 
its up state, the Husimi function i s sp read in the phase 
space after t = 20 kicks (left plot) [i^- On the contrary, 
if the qubit is in its down state, the Husimi function is 
localized around the stable fixed point at t = 20 (right 
plot). The difference between the Husimi distributions 
in these two cases is evident and leads to measurable ef- 
fects. Indeed, as shown in Fig.lHl the values of the second 
moment < > are very different at short times for up 
and down spins. 

It is interesting to discuss the case of weak system- 
detector coupling with e 1, ec = eft <C 1. In the chaotic 
regime for the kicked rotator the Husimi function exhibits 
hypersensitivity to perturbations 46J. Indeed, as shown 
in Fig. IIUI we have markedly different Husimi functions 
when the detector is coupled to up or down spin. In the 
case of Fig. UHl we have e = 0.4 and H = 1.23 x lO'^, 
corresponding to Ec = 4.91 x 10~^ ^ 1. However, it 
is unclear how to extract measurable information from 



FIG. 8: (Color online) Husimi function in action-angle vari- 
ables (p, 8) for the detector, with —n < p < n (vertical axis) 
and < 9 < 2t7 (horizontal axis) for the kicked rotator cou- 
pled to up spin (left) and down spin (right), at K — 4.5, 
Ec = 0.8, 5 = 0.1, h = 1.23 X 10"^, t = 20. The initial states 
of the kicked rotator and the qubit are a Gaussian packet cen- 
tered at the fixed point p = 0, 61 = tt and = (|0) + 1 1))/\/2. 
Color represents the density from blue/black (minimal value) 
to red/gray (maximal value). 
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FIG. 9: Time dependence of (p^) for up spin (curve) and 
down spin (dashed curve). Parameter values are the same as 
in Figil 



this difference. Indeed, it is reasonable to assume that 
only the coarse graining properties of the semi-classical 
detector are accessible, the size of the coarse graining 
being much larger than the Planck cell. The difficulty 
to distinguish the left and right Husimi plots of Fig. ^| 
after coarse graining is illustrated in Fig. ^2 lu this 
latter figure we compute the integral Wd of the Husimi 
distribution over a box of size 2.53 x 10^^ ^ h = 1.23 x 
10^^, centered, as the initial wave packet, at p = 0, 6* = 
TT. The evolution in time of Wd is shown for different 
but small coupling strengths {ec ^ I). It is clear that 
the evolutions of Wd (t) for up and down spin states can 
be hardly distinguished. 
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FIG. 10: (Color online) Hypersensitivity of the Husimi func- 
tion on the spin value for K — 8, e = 0.4, 5 = 0.2, and 
h = 1.23 X 10"^ From top to bottom t = 0, 4, 8, 12. The 
left plots are for up spin, the right ones for down spin. The 
initial states of the kicked rotator and of the qubit are a Gaus- 
sian packet centered at the fixed point p = 0, 9 = tt and 
\ips) = (|0) + \l))/V2. The color code is as in Fig.|Sl 



V. CONCLUSIONS 

In this paper we have proposed and studied a deter- 
ministic model of quasi-classical detector coupled to a 
single-qubit system. Our results show that the detector 
introduces dynamical decoherence in the system. More- 
over, we have shown that an efficient measurement is 
possible in the case of strong system-detector coupling. 
In the case of weak coupling, the chaotic dynamics of the 
detector is still hypersensitive to the state of the qubit 
but it is unclear how to detect the up/down states in a 
coarse graining measurement. It is an interesting ques- 
tion whether this conclusion remains valid when the un- 
avoidable coupling of the detector with a dissipative en- 
vironment is taken into account. 

We thank the IDRIS in Orsay and CalMiP in Toulouse 
for access to their supercomputers. This work was sup- 
ported in part by the project EDIQIP of the IST-FET 




FIG. 11: (Color online) Integral Wd of the Husimi distribu- 
tion over a square box of size 2.53 x 10~^ ^ h = 1.23 x 10~^, 
centered at p = 0, = tt, for spin up (dashed red/gray curve) 
and spin down (solid black curve) states. From top to bot- 
tom: (e = 0.4, 5 = 0.0), (e = 0.2, S = 0.2), and (e = 0.4, 
S = 0.2). The initial state of the detector and the other pa- 
rameter values are as in Fig. 1101 The top and middle curves 
are shifted in the j/-axis direction by +2 and -|-1, respectively. 



program of the EC. 



APPENDIX A: MASTER EQUATION 

In this Appendix, we provide an alternative derivation 
of the phase damping map (jl5(l . based on the master 
equation approach. For this purpose, we start from the 
overall Hamiltonian of Eq. Q and make the usual Born 
and Markov approximations. That is to say, we assume 
that the detector's dynamics is practically unaffected by 
the interaction with the qubit and that any effect the 
qubit has on the detector is limited to a time scale much 
shorter than the time scales of interest for the dynamics 
of the qubit. The first condition is fulfilled in the case of 
weak coupling, e ^ 1, the second is satisfied due to the 
fast decay of correlations for the detector in the chaotic 
regime (here we assume that 5^1 and that correlations 
decay in a few kicks). Under these assumptions, we can 
derive 0| the following master equation: 



dp 
'dt 



= -i[Hs,p] + y ((7z/5(72 - p) ^S{t - m). 



(Al) 



where we have used 



elTTd[cos''{e)pd 



E 



S{t — m) 



(A2) 



Here pd is the detector's density matrix and the tilde 
denotes the fact that we are using the interaction pic- 
ture (that is, the time evolution of the operators with a 
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tilde is ruled by the detector's Hamiltonian Hd)- Assum- 
ing that the density matrix corresponds to a uniform 
distribution in the 9 variable, we obtain Tidicos^ (9) pd] ~ 
^ J^^ d9 cos^ ^' = 5 ■ The integration of the master equa- 
tion HA1|I in one time step leads to the phase damping 
map l(T3|) . 

APPENDIX B: CONTINUOUS MODEL 

The continuous version of the phase damping map H15|) 

is 

X = — e^x, 

y = -e^y-2Sz, (Bl) 
z = 2Sy. 

The solution for x{t) is a simple decay with the rate 
r = e2: 

x{t) = a;(0)e"". (B2) 

The solution for y{t) and z{t) reads as follows: 

r y{t) ^ [asin(cjt) + bcos{Ljt)]e~'^*^^ , , . 

\ z{t) ^ [csinM) + dcos(tjt)]e-"/2^ ^ ' 

where 

^ = ^452-^ (B4) 



and the coefRcients a,b,c,d can be expressed as a func- 
tion of the initial conditions y(0),2;(0). Therefore, for 
e <C in the a^-basis both the diagonal and the off- 

diagonal elements of the density matrix decay as ^, in 
agreement with our numerical data shown in Sec. III. In- 
deed, we have at large times |pii — ^1 = 1^1 ^ exp(— rt/2) 
and Ipoil = + y2)i/2 ^ \y\ ^ exp(-rV2). 

When the couphng to the detector becomes strong, so 
that r > 4(5, then the oscillations in y{t) and z{t) turn 
into decay with two characteristic rates: 



r /f2 

r±= 2 ±y--4<52. (B5) 



The smallest rate (r_) describes the quantum Zeno ef- 
fect. For r ^ 45 we obtain r_ ~ AS'^/T, in agreement 
with our numerical results shown in Fig. |31 

It is important to remark that in the basis of the eigen- 
states of the Hamiltonian Hg one can recover the relax- 
ation and decoherence time scales usually discussed in 
the literature. In this basis x describes the deviation of 
the diagonal elements of the density matrix from 1/2, 
while y and z give the decay of the real and imaginary 
part of the off-diagonal matrix elements. Therefore, for 
e < the diagonal elements decay with rate F and the 
off-diagonal elements with rate F/2. 
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